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Models which allow an explicit application to structurally modulated substances are reviewed 
within the frame of a symmetry-based approach starting from discrete lattice theory. Focus is set 
on models formulated in terms of local variables assigned to discrete crystallographic units (unit 
cells or parts of them). Especially considered are symmetry-based pseudo spin models. Methods are 
discussed which permit the handling of the statistical mechanics of such models, the translation of 
model calculation results to a form allowing the prediction of experimental data, and the theoretical 
"^C," , determination of point defect influences. 
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I. INTRODUCTION 



An interesting class of materials exhibits, as a distinguishing feature, phase diagrams with several, in some cases 
with very many, commensurately and incommensurately modulated phases. At present there exists a large body of 
experimental data on a multitude of such materials, and various types of models have been proposed as bases for a 
^ ■ theoretical description. An assessment of their merits and their shortcomings depends i) on the weight which one 
attaches to the physical justification of the general form of the model and of the special values of its parameters, ii) 
on the properties of the real system which are to be explained, and iii) on the degree of agreement between theory 
and experiment which one requires. 

At the present state of research it seems appropriate to survey a large class of these models on the same footing, 
' to discuss how they are connected with the original atomistic or lattice theoretical basis, to find, if possible, a 
^\ , methodological and hierarchical classification scheme into which they can be fitted, and to compare model predictions 
of properties like phase diagrams, values of physical quantities etc. with each other and_with experimental data. 

There exist extensive discussions of models exhibiting spatially modulated structuresM. In the present paper it will 
be our special concern to discuss how different models are obtained when the basic first principles theory is gradually 
more and more simplified so that finally theoretical descriptions arise which allow explicit, though only approximate, 
calculations of physical properties. We will be particularly interested in various models whose common feature is 
that they describe the macroscopic specimens in terms of local variables assigned to discrete crystallographic units on 
q \ the atomistic level. We will be especially concerned with a class of such models, which arise if the continuous local 
variables are projected onto a number of discrete states and, finally, with a description of these states in terms of local 
pseudo spins with one or several components. 

Besides the derivation of model characteristics - e.g. the global structure of phase diagrams, detailed information 
on the formation of new phases by structure branchings and their accumulation points (which are connected with 
transitions from first to second order phase transitions) or the exact structure of phase diagrams in the vicinity of 
multiphase points - we shall also report methods for establishing a connection between theoretical (model) data and 
experimentally accessible properties of real systems. 

An interesting additional aspect of the theoretical treatment is the possibility to incorporate impurities into the 
theory. On the one hand, this line of the theory is of practical relevance for the understanding of real, imperfect 
crystals; on the other hand it is of interest to find out if the basic ideas of the modelling procedure carry so far as to 
allow also the description of impurity effects. 

The paper is organized as follows. After a brief reference to the experimental situation and some general remarks 
on the interrelation of model parameters and measured quantities in section [0], short comments on phenomenological 
theories of modulated structures and on the derivation of pressure-temperature or temperature-electric field phase 
diagrams from Landau type thermodynamic potentials are made in section |ni| In section pV| a symmetry-based lattice 
theoretical description of uniaxially modulated structures is presented which provides a framework for reviewing a 
numbe r of models in the subsequent sections (models with continuous variables in section [v[ pseudo spin models in 



section 



VI ) . Section VII describes methods for bridging the gap between the results of model calculations (in terms 
of model parameters) and experimentally observed properties (in terms of experimental control parameters) . These 
methods do not postulate an ad hoc correspondence but - with further approximations - make use of the lattice 
theoretical foundation of the models and of macroscopic material properties. The incorporation of point defects into 
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the theory is discussed in section VIII . In section IX , a summarizing discussion and an outlook to future developments 
is given. 



II. EXAMPLES OF MODULATED SYSTEMS; OBJECTIVES OF THEIR THEORETICAL 

DESCRIPTION 



During the late 70s and 80s a worldwide hunttra for materials with cascades of structurally modulated phases led 
to a large amount of experimental data on numerous compounds. The largest class of such crystals is formed by 
A'A"BX4-compounds, among which Rb^ZnCU is a prototypical member. Apart from measurements under different 
temperatures, pressures and electric fields, homologous series of A'A"BX4-compounds, e.g. the tetramethylammonium 
tetrahalogenometallate (TMA — MX4) series, allow the investigation of the influence of the chemical composition on 
the sequence of phase transitions. 

An outstanding specific material with an especially large number of commensurately or incommensurately modu- 
lated phases is betaine calcium chloride dihydrate (BCCD). Its phase diagram, which has been thoroughly investigated 
experimentally e.g. with respect to its dependence on temperature, pressure, uniaxial stress, electric fields, or impu- 
rities shows many very interesting features. The observed global and local details of the phase diagrams invite a 
comparison with theoretical predictions to test the performance of various models. 

We only brieJBy_mention three further special examples since surveys on structurally modulated materials were given 
elsewhere (e.g.ETu). These examples are interesting in that they show that in families of isomorpheous compounds 
materials with properties deviating significantly from the majority behavior may be observed: a) All investigated 
modulated substances of the TMA — MX4 family have their modulation wave vector q along a with the exception 
of [N(CH3)4]2CuBr4 with q parallel to c (similar to its homologuoJ tetramethylphosphonium-CuBr4). b) Rd^ZnCLi 
and Rb2ZnBr4 on the one hand and Cs2CdBr4 and Cs2HgBr4 on the other hand show cascades of modulated phases 
but Cs2ZnCl4 and Cs2ZnBr4 do not. c) Whereas most of the members of the A2BX4 family exhibit phase transitions 
of purely order-disorder type and hence lack a soft modeS, K 2 Se04 is primarily displacive in its incommensurate 
transition, but also possesses a significant order-disorder component. 

Most of the previously mentioned materials with uniaxial structural modulations were investigated at ambient 
pressure over large temperature intervals; for some substances, pressure was also varied or the influence of an electric 
field was studied. 

Different approaches offer themselves for a theoretical interpretation of these data: a rather rigorous treatment 
would start from the basic Hamiltonian formulated as a function of atomic positions, momenta, and externally imposed 
parameters (e.g. homogenous stresses and electric field components). This Hamiltonian might be deduced e.g. from an 
ab initio electron theory or from a formulation in terms of empirical interatomic or intermolecular interactions. This 
treatment would then employ the methods of statistical mechanics to derive for example thermodynamic potentials for 
equilibrium structures and phase diagrams. These would be obtained as functions of temperature and experimental 
control parameters and could be directly compared with experimental results. If (due to computational difficulties) the 
method is not (yet) feasible, one can make use of semiempirical treatments employing suitable models. In the standard 
procedure for such a model description of modulated systems, physical properties like phase diagrams are derived 
from effective model Hamiltonians with the help of approximate methods of statistical mechanics. The results are 
primarily obtained as functions of temperature and model parameters (including e.g. effective intercell interactions 
and electric field components). For a comparison with experimental data it would be desirable to find physically 
motivated connections between the effective interactions occuring in the model Hamiltonian and experimentally given 
parameters like temperature and applied stresses. 

Most applications of model calculations to the interpretations of empirical data restrict themselves, however, to a 
comparison of general topological features ftfj-the model phase diagrams with those obtained experimentally or try to 
establish an ad hoc correspondence. In ref£rE2l an explicit mapping from a theoretical temperature-interaction-phase 
diagram to a pressure-temperature phase diagram was presented which makes use of material properties like therma l 
expansion, elastic constants etc. The method was used exernplariU for BCCD and is described in subsection VILA. 

Phase diagrams were calculated. e.g. for the DIFFOUR modeMt4 the-ANNNI model and its extensionsMlia the DIS 



modeS^, the AANNDI models, and for Chen and Walker's modeltfll^. Usually two-dimensional planar parameter 
space sections were considered. 

All models discussed in the sections |y|-|v| contain variables that correspond explicitely to microscopic crystalline 
quantities: both local generalized coordinates or pseudo spins describe degrees of freedom of crystallographic (sub)cells. 
These model variables are coupled by effective interactions, which correspond to certain thermal averages of combina- 
tions of original atomistic interactions. If it were possible to find a temperature-, stress-, and electric field-dependent 
mapping of the atomistic couplings onto the model parameters, precise predictions for the occurence of phase tran- 
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sitions could be made. Particularly, it would be possible for a special substance to transform a theoretical phase 
diagram. 

Because of their experimental relevance, phase diagrams as functions of temperature and pressure are of special 
interest. The set of states experimentally accessible by variations of T and p corresponds to a curved two-dimensional 
surface in the higher-dimensional model parameter space spanned for example by temperature and effective interac- 
tions. In general, this surface will not be planar. 

III. PHENOMENOLOGICAL THEORIES 

As the phenomenological Landau theory has been the subject of various reviews (see for examplelli) we will focus 
in the following on different schemes for describing successive transitions between phases characterized by different 
wave numbers. We will especially discuss two approaches to handle a cascade of c©Hiiapnsurate phases. 

In order to describe the cascade of phase transitions of a devil's staircase SannikovE3c3 proposed a phenomenological 
theory in which the sequence of phase transitions is determined by one soft branch of the lattice vibrations. The 
thermodynamic potential for a transition to a commensurate phase characterized by the wave number qi = m/l can 
be written in the form 

$ = ap 2 + /3p 4 - a[p 21 cos 2l( (1) 

where p and £ are the amplitude and the phase of the normal coordinates. The coefficients a[ differ from zero only at 
the points of the Brillouin zone given by q = qi = y where m and I are integers, i.e. for every I one obtains a different 
potential <I> (the wave numbers q will be given in units of 2n/ai, i = 1, 2, 3, being the lattice constants). The limit 
I — ► oo corresponding to irrational values of q yields the potential 

$ = ap 2 + Pp 4 

describing the transition to an incommensurate phase. 

The coefficient a in eqn. (Jl|) is a continuous function of q. The approximation used for this function depends upon the 
interval of q- values-considered. If the end points of the interval do not satisfy the Lifshitz condition, a is approximated 
by the expressions 

a = a + 5 (q — b) 2 . 

Here b is the wave number of the soft mode responsible for the minimum of the branch a(q). This first case is realized 
for example in the TMA — MCI4 family. If, on the contrary, one of the end points satifies the Lifshitz condition (e.g. 
q = 0) the proposed form iso 

a = A + K(q 2 -B 2 ) 2 

which pjsults from an expansion in q 2 around the point q = 0. Examples for this case are thiourea [SC(NH2)2] and 
BCCEO. In the following we will focus on the latter case. 

Introducing new dimensionless variables x = ^ and y — -^s an d assuming that the value \a\\ is the same for all I 
Sannikov derived a phase diagram in the x-y-space. Q has the same dimension as q and its value is arbitrary. Under 
the assumption of a linear T-dependence of A and B (and therefore of x and y) the phase sequences were obtained 
for different materials. 

In a detailed analysis for BCCD, thermodynamic potentials for various values of q were derived. Using an ad hoc 
linear transformation from dimensionless model variables to T and p and choosing tentative values forptiie remaining 
parameters a rather satisfying fit of theoretical to experimental T-p phase diagrams could be obtainedHa. 

In order to study the influence of external (electric) fields E 1 , coupling terms between the fields and the two- 
component order parameter were added to the potential (^). The transformation properties of the conjugated quan- 
tities Pi = —-§^7 (spontaneous polarization) could than be derived and compared to experimented. 

A diffeijenL-phenomenological theory was developed especially for the interpretation of the phase transition sequence 
of BCCDE3 Lj. Noticing that the polarization of the ferroelectric low temperature phase (q = 0) corresponds to the 
primary order parameter, the starting point for this approach is the free energy density (in reduced units) 
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The order parameter is supposed to be site-dependent. The terms — j (§§J an< ^ \ y§^J describe the occurence 

of modulations whereas the term /i£ 2 ^§§^ is responsible for the temperature dependence of the wave vector0. In 

addition umklapp terms were introduced into the expression for the free energy density of the commensurate phases 
and approximated by a forth degree term in the order parameter amplitude, 

CAimklapp = -fait, 21 ~ --70e$(q) £ 4 , 

for a commensurate phase with wave number q = y. 

Using the plane- wave approximation (i.e. £ = £o cosfcx) and choosing appropriate values for /3qfh(q), the phase diagram 
was derived by computing the free energies of the different, phases as functions of temperatureEH. Introducing coupling 
termSr-between the order parameter and an electric fielcO the E — T phase diagram was calculated. In a subsequent 
paperEa the sinusoidal scheme was abandoned and the occurence of a multisoliton regime was considered. 



IV. DERIVATION OF MODELS FROM AN ATOMISTIC LATTICE THEORETICAL BASIS 



A. Formulation in terms of local modes 



Models based on the discrete structure of the lattice can provide a way to achieve a physical understanding of the 
processes on the atomic level of the crystal which goes beyond a mere phenomenological description. Such microscopic 
models can be handled analytically or numerically. In the subsequent sections models which have been closely studied 
with respect to the application to experimentally investigated substances are discussed. Each model variable - either 
continuous as in the case of the DIFFOUR models or discrete as in the case of pseudo spin models - usually describes 
one degree of freedom per lattice site. Previously it was considered a weak spot, especially of pseudo spin models, that 
the significance (on an atomistic scale) of the model (pseudo spin) variables was not well defined. For this reason we 
summarize a general procedure for the systematic definition of symmetry-based model variables on a lattice theoretical 
basis. This procedure stresses the atomistic roots of the models, leads to models conforming to lattice symmetry and 
facilitates the prediction-of spontaneous polarizations of modulated phases and of the .atomic displacements occuring 
at structural transitionsc2l. It combines the idea of local modes formulated e.g. in refEi with the rigorous symmetry 
considerations of Landau's theory (for instance, see refOEj). It leads to a Hamiltonian that serves as a starting 
point for various approximations leading to different models (e.g. the DIFFOUR models or the DIS model) and thus 
provides a frame for reviewing very different models from a common point of view. 

The key idea of the general procedure for the derivation of a microscopic model Hamiltonian of this procedure is the 
introduction of variables as generalized coordinates of symmetry-adapted local modes (S ALMs) . This facilitates the 
incorporation of both the discreteness of the lattice and the overall symmetry of the crystal. The symmetry properties 
of the SALMs can either be derived from first-principles calculations or taken from experimental data, i.e. from the 
observation of the symmetry of the displacements responsible for the phase transitions. 

Let the crystal be composed of N unit cells each containing K particles (of which K' are non equivalent) and 
assume, for definiteness, Born-von Karman periodic boundary conditions. The crystal structure is then described by 
a 3K iV-dimensional configuration vector R, where each entry specifies a component of the position of an atom. An 
arbitrary configuration 

E = Ro + SR 

is decomposed into a contribution R , which is invariant under all elements Q = {R|t} of the high temperature normal 
phase space group (S , and a symmetry- breaking displacement i5R, which is zero above the transition temperature 
T cr it from the normal phase to the lower symmetry phases (unmodulated or commensurately or incommensurately 
modulated). The displacement SK from a reference structure at given reference values of T, p, etc. can be written as 
a superposition 



5R = 5>R n (2) 



of local contributions <5R n . They describe displacements for only those atoms which are associated to unit n, having 
zero entries for all others. These units do not necessarily correspond to the unit cells of the crystal; in the case of the 



space group Pnma discussed in subsection [V C , half cells are chosen as units accounting for the pseudoperiodicity of 
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such crystals along one of the lattice vectors. <5R n is expanded in terms of local modes. In most cases it is sufficient 
to retain only a few (e.g. one or two) local modes per unit. 

It is useful to choose the local modes such that they reflect the crystal symmetries: the smallest unit in the crystal 
with respect to symmetry is the asymmetric unit 21q of the space group 25 - Symmetry imposes no restrictions on the 
position of the K' atoms associated to this polyhedral subcell. Further subcells 21^ defined by 

% = G 1 ^ and 2l l n = {E|n}% 

cover the whole crystal. The B space group operations Q 1 , . . . , Q B are coset representatives in a coset decomposition 
of ©o with respect to its invariant subgroup X of primitive translations. The index B — [C5 : T] equals the order of 
the point group ©o of O and determines the number of subcells 2l„ per cell n. The coset representatives have to be 
chosen such that these B subcells form a contiguous new choice of cell. Any displacement of the K' atoms associated 
with 21^ can be expanded in terms of a basis set of 3K ' subcell modes VJf , k — 1, . . . , 3K'. They should satisfy the 
relation VJ," = QV^ where Q € <3 is the space group operation transforming subcell 2l„ into 21^,. Thus, the modes 
in any subcell are entirely determined by the (arbitrary) choice of modes in the asymmetric unit 21q. 

A complete - for K'B > K over-complete (see below) - set of 3K'BN local modes can be constructed as linear 
combinations of the subcell modes 

B 

Wi K = ^YiT (3) 

i=l 

with a B x B transformation matrix S, the form of which is determined by the IRREPs of the point group (So- 
For fixed Ik each local mode W^ K describes the same displacement pattern of atoms in their respective unit, i.e. 
W^+n = { E l m }WjT- We will call the set {W^ : k = 1, . . . ,3K'} the 7th set of symmetry-adapted local modes 
for unit n. Any linear combination of displacement vectors in one set of SALMs has the same symmetry properties. 
Apart from a primitive translation, an arbitrary space group operation only permutes the subcell indices and thus 
the transformation behavior of the set of W^" is entirely determined by S. The actual positions of the atoms in the 
subcells have not to be known. Only those sets of SALMs can contribute to the structural deformation U that have 
the symmetry properties imposed by the IRREP according to which U transforms. 

The existence of atoms on special positions, i.e. on faces, edges and corners of the subcells causes the number 
3K'BN of local modes to be in general higher than the number 3KN of independent displacements. Hence, the 
coefficients of an expansion of a general displacement in terms of these SALMs have to satisfy consistency relations. 

The symmetry-breaking displacement (SR n in unit n is decomposed into contributions <5R„ from the B sets of 
SALMs: 

B B 3K' 

^ n = E^ = EE a n K In K - (4) 

7=1 J=l«=l 

Under a continuous change of external parameters the average positions of the atoms are displaced in such a way that 
a valley of low energy is followed in the high-dimensional configuration space. The symmetry-breaking displacement 
of any atom in the asymmetric unit 21q which is connected with the transition is not determined by space group 
symmetry. It will, in general, be different even for materials exhibiting the same symmetries of the normal and the 

modulated phases (e.g. A2BX4 compounds). Hence, the contribution 5R^ = Y^kLi n n K W^, K from the ith set of SALMs 
will, in general, lie on a d a < 3i^'-dimensional curved surface or, in its simplest and most frequent case, describe a 
curved path [d a = 1) upon changing temperature or pressure. Either case corresponds to a simultaneous variation 
of in general all coefficients a^ K in equation (|J) as the fixed subcell modes V„ can not simultaneously be adapted 
to the curved paths described by every SR^. Thus, instead of using the coefficients from equation (^) as generalized 
variables for the determination of the displacements in unit n, one introduces new local sets of curvilinear coordinates 
for each set of SALMs by means of a coordinate transformation a^ K = j lK {Q^ , ■ ■ ■ , Qn 3K ), « = 1, . . . , 3K' such 
that the contribution 

3K' 

SBL(Q?, • • • , Q^ K ') = E ■ • • ■ Q^ 3K 'Mi K 

K=l 

to the symmetry-breaking displacement in unit n [cf. equation (||)] occurs when Q^ 1 , . . . , Q{f a increase from zero to 
some non-zero value whereas all other generalized variables remain at their thermal average of zero. In the simplest 
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case d a — 1, only one relevant variable Q 1 ^ per set of SALMs is needed and thus the value of Q 1 ^ determines a 
linear combination of displacement vectors from the /th set of SALMs for unit n (which has the same transformation 
behavior as any other linear combination of SALMs from the same set). 

On the one hand, the space group symmetries and the direction of the modulation determine the subcell generators 
Q 1 ,. .. ,Q B and the matrix S and thus should be the same for a whole class of materials. On the other hand, the 
displacements of atoms in the asymmetric unit and thus the coordinate transformations f Ifi are determined by the 
local potential and thus by the specific material under investigation. This allows a separation of properties specific 
for a certain class of materials from properties specific for the substance under investigation. 

The vector field 8l£(Q™, . . . ) gives a contribution to the structural transformation in unit n with definite trans- 
formation behavior. In general, d s sets of SALMs (I = 1,... ,d s without loss of generality) are needed for the 
construction of symmetry modes transforming according to IRREPs of the group © (q) of the wave vector. When 
©o has only one-dimensional real IRREPs and only one order parameter transforming according to one IRREP of 
©(q) gives the primary contribution to the structural deformation (as it is the case in BCCD), only one (d s = 1) 
certain set of SALMs (I — 1) and the corresponding symmetry-breaking displacement (5R* is needed for every unit. 
The reverse is not true because by different spatial modulations of Q^ 1 one set of SALMs can be used to construct 
symmetry modes transforming according to different IRREPs of © (q) and different sets of SALMs may contribute to 
symmetry modes with same transformation behavior (e.g. at a point q ^ inside the Brillouin zone when ©(q) is a 
true subgroup of ©o; examples will be given in subsection IV C). Especially, the set of SALMs corresponding to the 
totally symmetric representation of ©o may be relevant because it may contribute to a structural deformation not 
transforming according to the totally symmetric representation of ©o- 

In the simplest case of one relevant set of SALMs (d s = 1) and one degree of freedom in the asymmetric unit 
(d a — I), expansion (0) can be written 



5R = ^<5Ri(Qi\0,... ,0) 



and the SifA-dimensional lattice theoretical problem is reduced to a system with one generalized local variable Q* 1 
per unit. In the next paragraph, for simplicity, we will restrict ourselves to this special case, which however can easily 
be generalized. 

The crystal potential is expanded in terms of the local generalized coordinates Q„ K . The terms containing only 
variables with given n form the local potential & loc (Ql^, . . . , Q B ' 3K ) for unit n. The general Hamiltonian of the 
system can then be written 



<& int (... ,Ql K ,... ,q£,...) 



where P^ K is a momentum canonically conjugate to Q^ K . Averaging out all irrelevant variables yields an effective 
Hamiltonian 



H = tr(pH) = ^2[(P^) 2 + t° c (Q^)] +r it (... ,Qi\...,Q 



' m 3 • 



(5) 



(The trace is taken only for the variables to be averaged out: I ^ 1,/t^ 1.) 

In the case of two relevant sets of SALMs per unit (d s — 2,1—1, 2), one gets the two-variable version 



n = 



£ [( P n) 



—hit 



(...,Q L n L ,Q 



11 ^21 
n ' 



o 11 o 21 



The same form with the superscript (21) replaced by (12) is obtained in the case of one relevant set of SALMs per 
unit (d s = 1, I — 1) but two relevant generalized coordinates (d a = 2) per asymmetric subcell. We will omit the 
superscript referring to n — 1 for brevity in the case d a = 1 of one generalized variable for each set of SALMs. 



B. Relation to other models 



The models discussed in the subsequent sections can be viewed as special cases of the Hamiltonian (j|) or its 
two- variable versions. The total number of relevant variables per unit is d a ■ d s . The local potential $ and the 
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interaction potential $ are effective potentials which, by their very definition, depend directly on elastic stresses 
(which determine the lattice constants) and on temperature due to the thermal averaging over the irrelevant variables. 
The potentials are usually expanded up to forth or second order respectively and only a subset of interactions between 
nearest, next nearest and third next nearest neighbors is kept. 

As in Landau theory, space group operations can be considered to act on the generalized variables instead on the 
modes W^. Hence, symmetry considerations give conditions for the dependence of the potential on the generalized 
variables, which have for example been taken into account when formulating the AANNDI or the DIS model. Moreover, 
the local potential may have equivalent minima: assume that there is one global minimum of the local potential for a 
certain configurational set of generalized local coordinates; then every space group operation which transforms a unit 
into itself but changes the configuration of the respective coordinates may generate another equivalent minimum. 

A further simplification may be introduced by projecting the remaining continuous generalized local variables onto 
two- valued pseudo spin variables as the signs of the former. This leads to different pseudo spin models: we discussed 
either the ANNMjnodel in the case of one relevant mode per half cell or the DIS model in the case of two relevant 
variables earliercHlE3. 

The Hamiltonian (^|) corresponds to a version of the DIFFOUR model with one degree of freedom per lattice site or 
to the ANNNI model. The more general case where two local modes with different transformation behavior (d s = 2) 
are incorporated corresponds to Chen and Walker's model or to the DIS model. Assuming two relevant displacement 
vector fields for the asymmetric unit (d a — 2), i.e. two generalized variables Qn 1 an d Ql? for the relevant set of 
SALMs, one arrives at a modeLcarresponding to the version of the DIFFOUR model with two degrees of freedom 
per lattice site discussed in ref-E-H (Q^ 1 and correspond to variables describing the displacements of the betaine 
and calcium chloride groups of BCCD respectively). For the description of BCCD with space group (So = Pnma, the 
units correspond to half cells. The latter two cases yield Hamiltonians with two different variables for each of such 
units. 

It is noteworthy that the Hamiltonian (Q) contains only one symmetry-breaking contribution per unit. As discussed 
above, the totally symmetric IRREP of <5q, which also contributes to the deviation of the displacements from the 
structure of the para phase at the original reference values of T, p, etc., is always present but has been absorbed into 
the atom positions of a new reference state. The latter can be viewed as giving a background contribution to the 
structural deformation of the crystal which is present throughout the whole range of temperature and pressure. 



C. Exemplary study of the procedure for materials with space group Pnma 



Among the 66 crystals with structurally incommensurate phases listed by CumminsS, the 22 materials belonging to 
the space group (So = D^h (in their normal phase) form the largest class. The materials in the A2BX4 (K^SeC^) family 
(apart from a small group of tungstates and molybdates), betaine calcium chloride dihydrate (BCCD), ammonium 
oxalate hemihydrate (AHO) and thiourea all exhibit the symmetries described by that space group. Depending on 
the particular choice of axes, various Hermann-Mauguin symbols have been attributed to these substances. For a 
more detailed discussion see ref.H. In the present study, we choose the assignment Pnma und restrict ourselves to 
modulations along c. Special attention is payed to the pseudoperiodicity. The reader should not get confused by 
the formalism as in the present case units correspond to half cells but the labeling with cell subscripts n is kept for 
simplicity. 

The asymmetric unit %l = [0; |] x [0; \] x [0; 1] of the space group Pnma is a parallelepiped (, block') containing 
all K' nonequivalent atoms from which the whole crystal can be filled by the B — 8 space group operations 

^{eiooo}, e 2 -K|oio}, g 3 = {C||lil}, g 4 = {i|m} 

and primitive translations {Ejniei + n^&i + 71363}. The eight blocks (see figure generated from the asymmetric 
unit by the operations G 1 , - ■ ■ ,G 8 are considered to form a new choice of cell at n = (000). An arbitrary block is 
labeled by the index i of the space group operation and the vector n of the primitive translation it was generated by. 
Blocks 1, . . . ,4 extend from z = to 773 + 1 (first half of cell), blocks 5, . . . ,8 from z = rz.3 + 1/2 to + 3/2 (second 
half of cell) along the c-axis. 

Local modes W„ K are constructed from the subcell displacement vectors V.n by means of the 8 x 8-matrix (cf. eqn. 
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Although constructed for the whole (new) unit cell, these local modes describe displacements for either the first 
(/ = 1, 2, 3, 4) or the second half (/ = 5, 6, 7, 8) of the cell when viewed along the c-direction in accordance with the 
definition of the blocks. Thus half cells (quadruples of blocks) are the units used to introduce the local potential. 
Nevertheless, we do not relabel the half cells but distinguish between front and back half cells by means of the subscript 
n and the superscript /. 

In Landau theory, i5R is decomposed into symmetry modes e kaA of the crystal transforming according to irreducible 
representations (IRREPs) T ka of <S : 



dim T K 



E EE E 



c kQA e ko 



(6) 



star k a 



where the first sum runs over all distinct stars, a labels the different IRREPs connected with star k, A accounts for 
the fact that the IRREP T kQ can appear more than once and the sum over s extends over all basis function of the 
IRREP r kQ . Contributions from the totally symmetric IRREP of <3q are included in R and excluded in equation 
(^) by the definition of SR. Below T cr i t at least one IRREP gives a non-zero contribution to equation (|TJ) and the 
respective symmetry-breaking displacement determines the structure of the newly formed phase. In the microscopic 
theory presented in this section, the symmetry modes c kaA e kQA are formed by superposition of the local mode vector 
fields <5R n . The generalized variables Q„, . . . are spatially modulated order parameters. 

It follows from the transformation behavior of the local modes under arbitrary elements of the space group ©o that 
it is possible to construct symmetry modes transforming according to the IRREPs A2 or A3 (but not Ai or A4) if one 
takes into account local modes with superscripts I = 3, 4, 7, 8 and transforming according to A+iOr A4 (but not A2 or 
A 3 ) if one takes into account I = 1,2, 5, 6 (0 < q < i, q = qc* , labeling of IRREPs as in refrj). For BCCD, which 
exhibits a soft phonon branch of A3-symmetry, the primary contribution is produced by the SALMs {W^, W„ K }. 
Additional secondary contributions (transforming according to the IRREPs A3 and A2 which were observed in the 
fourfold phaseEH) are accounted for by a second set of SALMs for each unit (half cell) ({W* K } and {W„ K } respectively). 
For the description of the relevant main contribution to structural deformations in BCCD, it is sufficient in both cases 
to keep one generalized variable for each set of SALMs, i.e. one has either two (Q^, Q„) or four generalized variables 

(Qn: Qli an d Qai Qn) P er ceu - 



D. The introduction of the pseudo spin formalism 

The effective potential <!> is a sum of local half cell potentials $ loC and harmonic two- variable couplings. To describe 

a system with 2 r equivalent minima of $ ° c , one keeps r relevant generalized coordinates Q^, . . . , Q r n per unit. The 
positions of the equivalent minima differ only in the signs of the r corresponding variables. In general, there are B 
equivalent minima per unit at the most. The ratio of the coupling strength and the height of the local potential 
barrier determines whether the system is more of the order-disorder or more of the displacive type. 

A pseudo spin formulation is derived in the following way: The interactions are expanded up to second order terms 
in the variables 

Ql = |Q£|sign Ql =: |Q>£. (7) 

These are replaced by the averages |Qn|c£ over single wells (taken separately for each |<2£|). Here 

cr£ := sign Q p n 

is just a shorthand notation for the sign of The interaction between the variables and becomes 

Kl ■ W\ ■ JQU ■ « =■■ JSnOga^, 
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with pseudo spin interactions = $ nm ■ \Qa\ ■ \Qm\- Since they are defined as thermal averages the latter depend 
quite naturally on temperature and stresses. The procedure yields the following pseudo spin Hamiltonian: 

W =^EE J n»m- (8) 
nm pq 

A switch of <t£ from +1 to —1 or vice versa is connected with a collective motion of all particles in the unit according 
to the corresponding SALM W„ which is associated with a transition from one well to another. Because a assumes 
only the values +1 or —1, it is called pseudo spin. It should not be confused with usual spins nor is it localized at a 
definite position in the unit—An association of pseudo spins with local normal coordinates for the description of soft 
modes was described in ref.ES. 

In this way, the distribution of the variables Qn over the wells in the unit is described by of the operators <r£. 
The formulation takes the actual average displacements and the thermal fluctuations in .the .neighborhood of a given 
equilibrium position into account. A more detailed analysis extending previous treatmentscJO to uniaxially modulated 
materials and incorporating both order-disorder and displacive phenomena will be reported in a futuiie publication. 

It should be mentioned that part of the SALM model was inspired by an idea presented in a papercil describing the 
distortions of four crystallographically equivalent ammonium tetrahedra in ammonium hydrogen oxalate hemihydrate 
by four pseudo spins. This model exhibited the possibility of incorporating symmetry considerations into pseudo spin 
Hamiltonians, although the interactions were not properly defined. 



V. MODELS WITH CONTINUOUS VARIABLES 



Many models proposed so far by various authors for the description of materials with modulated structures can be 
connected with the atomistic basis along the lines discussed in section |v|. In the present section we shall refer to 
some models with continuous local variables. 



A. The Frenkel-Kontorova and DIFFOUR models 

In the Frenkel-Kontorova modelilS a linear chain of atoms with harmonic nearest neighbor couplings in an external 
periodic potential incommensurate with the underlying lattice is considered. The incommensurability of the external 
potential (period b) with the lattice distance a in the potential 



V = f E - ^ - a) 2 + * ' (1 - cos(2tt^/6))] 



of the Frenkel-Kontorova model produces the frustrations necessary for the appearance .of (in)commensurate struc- 
tures. In the continuum limit of this model (also known as Frank-van der Merwe modelo) the equations of motion 
lead to the sine-Gordon equation of which the solutions correspond-to a soliton lattice. 

In the beginning of the 80s, Janssen and Tjon proposed a modelc^a that they later called , discrete <I> 4 ' (DIFFOUR) 
model following a similar model in field theory. Ref.O provides an extensive discussion. The DIFFOUR model is 
intended as a simple one-dimensional model featuring most of the typical properties found in modulated crystal phases. 
Its essential ingredient is the competition between short-range interactions. In its simplest version it introduces a 
linear chain with anharmonic nearest neighbor and harmonic next nearest and third nearest neighbor interactions. 
The frustration mechanism is the competition of the next and third nearest neighbor interactions. The reduced 
potential of the translation invariant version is given by 



^rcd = 



a . . 2 (3 . . 2 8 . 

— (vn-ut-i) +-{ut-ue-2) +2*- U£ 



u t - 3 ) 2 + ~{ut - ut_i) 4 



where ug denotes the displacement of the Ith. particle in the chain, ug is scaled such that V^ e d contains only the 
necessary model parameters a,f3 and S. Substituting xg := ug — ug-\ yields the not-translation invariant version: 



V Ied = 



A 



1 



Bxixi-i + DxgXg-2 



This is the potential of an effective Hamiltonian with a local 2-4-potential and nearest and next nearest neigbor 
interactions which can be considered as the result of thermal averaging over the non-relevant variables, while the 







quantity xg is interpreted as the relevant degree of freedom for position £ in the chain. The coefficients the depend 
quite naturally on temperature and other control parameters. The stable structures fulfilling dV rc d/dxi = and 
having the smallest free energies represent the ground state configurations. A temperature can be introduced either 
by identifying the potentials as the respective thermodynamic potentials (e.g. free energies) and assuming that the 
parameters of the model are temperature-dependentc3 or by starting from a Hamiltopian with the above potential 
and treating (a three-dimensional version of) the model in mean fieid approximatiorO. 

An extension of this model to twaxlegrees of freedom per positior£3 leads to a DIFFOUR model for the description 
of the phonon spectrum of BCCDoo. 



B. Chen and Walker's model 



Contrary to the approach of the DIFFOUR models starting from a general formulation and adapting the coefficients 
to specific substances later, Chen and Walker use information about the crystal structure and the symmetries of 
excitations in uniaxially modulated materials whose high temperature phase exhibits the symmetries of the s paic e 
group D\\ from the beginning. Their model was adapted to the class of A 2 BX 4 compoundgl3'E3 and to BCCDc§E3. 
The crystal is decomposed in equidistant layers with a layer spacing of half the lattice constant along the direction of 
modulation. Displacements of the atoms in layer £ are given byE2l 

veei(r 2 ) + w e ei(T 3 ), 

where ei(T2) and ^3) are layer modes belonging to the irreducible representations T2 and T3 of the two-dimensional 
group Civ , respectively, from which three-dimensional symmetry modes transforming like the observed soft mode can 
be constructed. Neglecting other layer modes, a Landau-type reduced free energy 



F = 



1 2 



1 4 

~A V l 

4 1 



1 / 2 

2 a H 



1 4 



~^2(Jvev e+1 + J' 



bviwj 



2^ 



(vewi + i ~ v t+1 wi) 



invariant under the transformations of the normal phase space group F>W with temperature-dependent parameters 
is introduced. Besides a local energy term, only nearest neighbor layer interactions are considered. The competition 
between the interactions J and J' on the one hand and the antisymmetric interaction term vgwz+i — V£ + \W£ on the 
other hand leads to frustrations and hence to the occurrence of commensurately and incommensurately modulated 
structures. 

Phase diagrams are calculated as planar two-dimensional sections through the higher-dimensional [a, a' \b, J, J')- 
parameter space: the specific structure of the layer variables {v£,wi} minimizing F determines the stable phase at 
every point of the parameter space. The respective space group can be determined from a given profile of amplitudes 
making use of the definition of ee(T2) and eg ^3). Up to three different space groups can be found for a given wave 
number of the modulation. A path is drawn in the parameter space of the theoretical phase diagram in such aojray as 
to obtain the same sequence of phases as determined experimentally as a function of T (for p = Q) for BCCBca. The 
expressions for the spontaneous polarization as expansion in terms of the layer variables can be used to determine the 
sequence of phases under an external electric fieldo. Observing that the soft mode in BCCD is an acoustic moder-a 
further layer mode coupled to the two optic modes and belonging to the irreducible representation T3 is introduced^ 
in the free energy. Phonon branches calculated from the latter show a softening of the lowest-lying acoustic mode for 
an appropriate choice of the eight reduced parameters. 



DIFFOUR mode l and ANNNI model (subsection VIA) on the one hand and Chen and Walker's model and DIS 
model (subsection VI C| ) on the other hand are closely related since in each case the pseudo spin model can be 
considered as the Ising limit of the model with continuous variables. 



VI. PSEUDO SPIN MODELS 



A. The ANNNI model and its extensions 
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1. The model 



The Axial Next Nearest Neighbor Ising (ANNNI) modell! was the subject of different excellent review articles00. 
In the following we will focus on the three-dimensional version and present the main results as well as some techniques 
used to examine the properties of this model. In the ANNNI model only one pseudo spin per unit (ijk) is considered. 
The Hamiltonian (without an external field) is given by 

T~L = —Jp/^^ijk + &i(j+l)k) — Jl <Tijk&ij(k+l) ~ J% / ] Vijk&ij{k+2} ■ (9) 

ij k ij k ij k 

Jo > is the ferro-type in-plane interaction between nearest neighbors whereas J\ and J 2 are the couplings between 
nearest and next nearest neighbors in direction of the modulations respectively. Frustration effects and therefore 
modulations arise for J 2 < and either sign of J\ . Due to the invariance of the Hamiltonian (^) under a simultaneous 
transformation J\ — > —J\ and er^ — > (— 1) CTyfe, it is only necessary to examine the case J\ > 0. Using reduced units, 
the ground state is given by the ferro phase for n < | and the (2}-phase for k > \ with n = — For the (2)-phase 

the spin structure repeats itself after four layers and is given by + H , i.e. two layers where all the pseudo spins 

have the value +1 followed by two layers with —1 pseudo spins. The two ground phases are separated by a multiphase 
point (MP) where an infinity of phases is degenerate. ,— ■ ,— ■ 

The A3NNI (Axial Third Nearest Neighbor Ising) modeE3~Eil differs from the usual ANNNI model by considering 
a third neighbor interaction in c-direction. This extension leads to a two-dimensional ground state phase diagram 
where some of the phases are separated by multiphase lines (ML). The two phases (12) and (3) occupy now a finite 
stability region at T = 0. r~in 

Starting point for the ELII (Effectively Long-range Interaction Ising) modelE3'E3 is the coupling between two degrees of 
freedom (one degree of freedom is considered to be an Ising spin) which leads to indirect effective long-range oscillating 
interactions J r = Y1 Jq eiqr between the pseudo spins, r being the number of layers separating the interacting pseudo 

spins. 



2. The mean field approximation 

In the mean field approximation (MFA)0E1 it is assumed that every spin interacts only with a mean field produced 
by the thermal averages of all spins interacting with it. This is equivalent to neglecting spin fluctuations or approx- 
imating the density operator by a product of single spin density operators. The free energy of the ANNNI model 
derived in this approximation (apart from constant factors) 

F = ^ [— 2J Sl — JiSkSk+i — J2SkSk+2 

k 

+k B T ((1 + S k ) In (1 + S k ) + (1 - S k ) In (1 - S k ))} 

is a function of the the mean field value Sk — {&ijk) of the fc-th layer. Minimizing the free energy leads to the 
equilibrium equations 



S k = tanh [4J S k + Ji (5 fc _ x + S k+1 ) + J 2 (S k - 2 + S k+2 )] . (10) 



The critical line 



fcsT , cr it(t) 



4J + 2(1 - k) Ji k<\ 
K 4J + (2/t+i)Ji k>\ 
which separates the disordered para phase from the ordered phases and the critical wave vector 



[0 K< \ 

<7crit = \ (11) 

[ arccos ^ k > j 

are obtained by linearizing the equilibrium equations ([l0|). At the Lifshitz point {n = \) the para phase, the ferro 
phase and the modulated phases meet. Figure (0) (from rcf.Ea) shows a typical mean field phase diagram with the 
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main commensurate phases derived for the case Jq = J±. Mean field phase diagrams of the A3NNI0 and of the ELII 
modeE3 have also beerupublished. 

Duxbury and SelkeE3'E3 verified for the ANNNI model the existence of structure combination branching processes 
responsible for the creation of new phases at finite temperatures. In this process two neighboring stable phases (A) 
and (B) produce a new phase (AB) at a definite branching point. At higher temperatures this process repeats itself 
leading for example to the appearance of the phases (A n B) at different branching points. The order of the phase 
transition changes at the accumulation points of the branching processes from first to second-order. The coordinates 
of the accumulation point of the (23") phases were determined by numerical extrapolations Using a fixed point 
expansion_thc. positions of accumulation points were calculated numerically and analytically for very different series 
of phasesET lj. This fixed point expansion will be the subject|p£]the next section.—. 

MFA seemed to indicate the existence of partially disorderecSO and asymm^icEZI phases for weak in- layer coupling 



Jo (i.e. Jo < 0.3Ji). In order to further examine this behavior NakanishitS used the mean field transfer matrix 
method (MFTM) where the interactions in c-direction are treated exactly, while the interactions in the layers are 
approximated by mean fields. The problem is reduced to one dimension with the free energy 



JqSI 



ksThiZ 



where the partition function Z = Tr exp 



can be calculated by use of transfer 



Ik J 

matrices, is the mean field of the fc-th layer. Minimizing the free energy with respect to the calculated thermal 
average S k leads to a system of nonlinear equations which have to be solved selfconsistently. As a result Nakanishi 
showed that the free energies of the partially disordered and the asymmetric states are always higher than the free 
energy of the symmetric state and concluded that the former states were an artifact of the MFA. Monte Carlo 
calculations^] supported this claim. 



3. The fixed point expansion 



The fixed point expansion method was developed for the calculatimi »f profiles, free energies and interaction energies 
of discommensurations in the frame of the ANNNI modetlfQ'u'LJO. The mean field equilibrium equations ([lOj) are 
reformulated as a four-dimensional mapping 

pfe+i = h(p k ) for the vectors P k = (S*_ 2 , S k - 1} S k , S k+ ijB. The pseudo 
spin profile P 7fc of a phase 7 with the periodicity N verifies the relation P 7fc = p^ k + N = h N (P_^ k ) and corresponds, 
thus, to an iterated fixed point of the mapping. 

Since a discommensuration is a localized defect the changes in the pseudo spin profile are the largest in the immediate 
vicinity of this defect. This region is treated numerically by means of the original (not linearized) equilibrium 
equations. The large number of pseudo spins whose numerical treatment would be cumbersome is conveniently 
handled by an expansion in the asymptotic region, however, about the fixed point leading to the linearized mapping 
A fe+1 = H fc A fc for the deviations A* = P 7fc — P fc . The matrices H fc are symplectic matrices, whose special properties 
allow a largely analytic derivation of the profiles, the free energies and the interactions of such defects. Of special 
importance is the (symplectic) product G k = Y[™ =1 ll k+n ~ tl where the period n of the matrices H fe (i.e. B_ k+n = H fe ) 

is n = N or n = y depending on the symmetry of the profile P 7fc . 

The asymptotic behavior of the deviations is determined in terms of the product matrix G k in the form A k+n = G k A k 

and therefore by the eigenvalues \ v of . Figure (||) shows the different fixed point types resulting from the positions 
of the four eigenvalues in the complex plane. The approach to the fixed point is thus monotonic for fixed point type 
2 (2') and oscillatory for type 3 (see ref.Ell for a detailed discussion of the fixed point types). 

Analytical results for the free energy S of one discommensuration or the interaction energy W2(d) of two discom- 
mensurations (which depends on the distance d between the two discommensurations) can be derived. If the fixed 
point type is 2 (2') the interaction W2(d) Letween two discommensurations is a monotonic function of d. For the fixed 
point type 3 the interaction is oscillatoryEl 

The fixed point expansion is very useful for the determination_of ±he order of the transitions and for a (numerical and 
analytical) calculation of the positions of accumulation pointsElTLj. Suppose two neighboring stable commensurate 
phases 7 and 7' are separated by a transition line. The phase 7' can be viewed as resulting from a periodic arrangement 
of Z adequate discommensurations in the phase 7 separated by d\,. .. ,dz (d^+z = d u ) layers. The difference of the 
free energies per layer is 
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z 

F = F i _ F i' = ZE + W(d u ... ,d z ) ZZ + J2w 2 (d v ) 

v=\ 

where we use the fact that usually, due to the rapid exponential decay of Wi (d) , only nearest neighbor discommensu- 
ration interactions are to be considered. The phase transition between the phases 7 and 7' takes place when F is equal 
to zero. If, at the transition, the interaction is monotonic [fixed point type 2 (2')] the transition is of second order 
and occurs where E is zero, the distance d of neighboring spontaneously formed discommensurations being infinite. 
If, on the contrary, the interaction is oscillatory (fixed point type 3) the transition is of first order and takes place for 
positive S. The distance of neighboring discommensurations is then finite. 

For a phase 7 three different lines meet at the accumulation point: the phase boundary, the line £ = and the 
line separating the fixed point types 3|-and 2 (211. This enabled Tentrup, Jenal, and Siems to find numerically the 
accumulation points of the series ( 9 , fc 3)H (23 fe )B and (2 fc+1 32 fc 3)Ei A low temperature mean field analysis allowed 
to derive analytical approximationsEa^o which render the coordinates of the accumulation points well. 



4- The low temperature series expansion 



In their seminal work Fisher and SelkeL30 derived the low temperature behavior of the ANNNI model by ana- 
lyzing an exact series expansion of the free energy. Herefore, they used the concept of structure variables l^, thus 
characterizing each possible ground state by an unique set of these variables. The variable isxlfifined by = 
where is the number of band sequences of type fi in a given structure in a lattice of L layeralS. The free energy 
can then be written in the form 



F = /2 



where the sum is over all possible band squences. The parameters a a are expanded in the form 



(k) 

the contributions resulting from the fcth order term in the series expansion of the free energy. In the fcth order 
term of the free energy enter all terms where the values of k spins differ from their gcoiiad state value. 
In their analysis of the different orders of the free energy Fisher and Selke concludedlZ§E3 that the infinity of phases 
with the phase symbol ^2 fe 3) are stable at low temperatures. Later Fisher and Szpilka refined this expansion using a 
transfer-matrix methocQ. They showed that among the periodic phases (2 fe 3) only those with 1 < k < fc max appear 
at any fixed small temperature. Since k max diverges as T — ► all phases of the form (2 fe 3) indeed spring from the 
MP. Similarly they concluded that the mixed phases (2 k 32 k+1 3) appear for feW(T) < k < fc( 2 '(T). Since — > 
and ► 1 in the limit T — > phases (2 fc 32 fc+1 3) appear arbitrarily close to the MP. 

This low temperature series expansion technique has also been applied to other models exhibiting modulations—for 
example to the A3NNI modelEJ, the three-state chiral clock ppdelEu, the generalized p-state chiral cl©ek modellI3, a 
six-state clock model with next nearest neigh bor interactional, or the DIS (Double Ising Spin) modelO. This latter 



model will be the subject of subsection VI C . 



5. The high temperature series expansion 
In the high temperature series expansior0~il the wavevector dependent susceptibility 

3 

is considered, where the sum is over all correlations and Zj is the number of lattice spacings in c-direction. Writing 
the correlation function 

<<W - Tr(e-^) 
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as an high temperature expansion one can construct a ^series for any wavenumber. As a result the Lifshitz point 
(for the case Jo = Ji) was found to occur fcpyi = 97P 2 H S3 I The critical wavenumber and the critical temperature 
as function of k were also calculated. In ref£-2 the critical exponents for the different regions along the para phase 
boundary (i.e., for the para phase- ferro phase transition, the para phase-modulated region transition, as well as for 
the multicritical Lifshitz point) were calculated for the ANNNI model formulated on_a.sc or fee lattice using a series 
expansion. It was found, in agreement with the results of renormalization group theoryai___, that the critical exponents 
of the phase transition from the para phase to the modulated region are given by the exponents of the 3d-XY model, 
whereas the critical behavior of the transition from the para to the ferro phase is that of the 3d-Ising model. 



B. The AANNDI model 

The Axial Antisymmetric Nearest Neighbor Double Ising (AANNDI) modelEl has originally been formulatedEl as 
an order-disorder mechanism of phase transitions for the compounds of the A'A"BX4-family. Two pseudo spins were 
introduced, describing the orientational states of a BX4 tetrahedron. The four discrete orientational states are labeled 
by o = ±1 (apex up or down) and t a = ±1 (turn to the right or left). Remarking that the pseudospins t a , a A 
and the product t a <j a transform according to three different onedimensional representations of the local symmetry 
group, the general Hamiltonian 

n-f T' \ "* „AA , tsI \ " _A_A 1 T i \ A A A A , T \ " „AA 

nn in nn in nn in nn out 

1 (7 V n-A^A , t ^AAAA , 1 «w \ " „AA I A A\ 

nn out nn out nn in 

+ l -M rtrf^-a,) (12) 

nn out 

was derived by retaining only symmetry-allowed nearest neighbor couplings between t a , a A and t a <j a . The sums 
nn in and nn out are over nearest neighbors in the planes and perpendicular to the planes respectively. 
The AANNDI model is derived from eqn. (jl^) by retaining the terms 

h = j$>V+i +^E^+! + ^E r ^+i - °t+0 (13) 

I I I 

where / labels the different layers. In-plane couplings assuring ferro or antiferro orderings in the planes are assumed 
but not taken explicitly into |account. The Hamiltonian ( |l3| ) ^therefore only a onedimensional one which is examined 
in mean-field approximationEj. In this analysis the producto u> A — T A a A was taken to be independent of t a and 
af thus leading to a mean-field free energy depending on the three layer order parameters (t a ), (rf) and (u A ). A 
closer look reveals that in the equilibrium equations (erf) decouples from the other order parameters. Therefore the 
term proportional to J can be neglectedlifl. 

It should be noticed that, retaining the antisymmetrical interactions in the planes (cf. eqn. (|l2|)), one obtains the 
PANNDI (Planar Antisymmetrical Nearest Neighbor Ising) model suitable for the investigation of twodimensional 
modulations in A'A"BX4-compounds. 



C. The Double Ising Spin (DIS) model 

The Double Ising Spin or DIS model has been proposed by u as a simple model for the investigation of 

uniaxially modulated materials with the para phase symmetry group D^- Its Hamiltonian can be derived from the 
general Hamiltonian ( |l2"| ) by introducing the two spins r = t a and a = <j a t a . Retaining only symmetric in-layer 
couplings (as we are not interested in twodimensional modulations) we finally obtain the DIS Hamiltonian 

ft = K'Y T ijkTij(k+l) + L'Y' J ijk' J ij(k+l) + ~2 E ( <7 *J'* T *J(*+1) ~ T i3k a ij{k+1)) 
ijk ijk ijk 

+ J E Tijk (f(i+l)jfc + Ti(j+l)k) + J' E ai 3 k ( a (i+l)3k + a i(j+l)k) (14) 
ijk ijk 
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Comparing it to the AANNDI Hamiltonian ( |l3| ) reveals that the unimportant term J has been dropped whereas the 
term L is considered. Furthermore in-layer couplings are explicitly taken into consideration giving way to a three- 
dimensional model. In the mean- field treatment we have only two order parameters per layer ((r) and (a)) compared 
to the three order parameters of the AANNDI model. 

Frustration and therefore modulations arise in c-direction because of the antagonistic effects of the symmetric nearest 
neighbor interactions K and L on the one hand and the antisymmetric interaction M on the other hand. The in-layer 
couplings J and J' are ferro couplings. 

In the ground state phase diagram (figure ^J) five different phases with a twodimensional stability region are separated 
by multiphase lines at which infinitely many different phases are degenerate. In the following we will only consider 
the multiphase lines K — L = — M (line 1) and K + L = —M (line 2). The lines 3 and 4 are obtained from the lines 

1 and 2 by the symmetry of the Hamiltonian. 

The global temperature-interaction phase diagram is examined in the mean-field approximation. One obtains very 
different twodimensional sections of the phase diagram depending on whether, at T — 0, the line 1 (figure^) or the line 

2 (figure ^) is intersected. In the first case the modulated phases are only encountered at rather high temperatures, in 
the second case all examined phases reach to low temperatures and seem to converge to the multiphase point. Only 
some of the commensurate phases are shown, the high-commensurate and the incommensurate phases filling up the 
shaded areas. 

Using .the fixed point expansion (cf. subsection VI A | ) we examined the processes leading to the appearance of new 
phasealj. In the case 1 the existence of structure combination branching processes and of accumulation points of 
these branchings has been verified, i.e. new phases are created by the same processes as in the ANNNI model. In the 
case 2 no branching points nor accumulation points have been found. The phase creation process is therefore different 
from the process in the ANNNI model. |— ■ 
In order to further clarify this different bphavior we performed a low temperature series expansion of the free energytZI 
following the ideas by Fisher and SelkeHj'LJ for the investigation of the ANNNI model (i.e. of a model with one- 
component spins) . This expansion reveals that in case 1 only four of the infinity of phases degenerate at the MP are 
stable at low temperatures. All other modulated phases are created by branching processes at higher temperatures. 
For case 2 it follows from the low temperature expansion that an infinity of different phases spring from the MP. 
At higher (but still small) temperatures some of the high-commensurate phases disappear. A close inspection of the 
matrices and vectors entering the computation of the free energiesEJ reveals that in the limit T — > the DIS model 
exhibits a complete devil's staircase for case 2. 

The critical properties of the model can be derived by writing the Hamiltonian (|l4|) in matrix form. A special 
local transformation thea leads to a Hamiltonian where the interaction terms are given by scalar products between 
two-component vectorsEJ. For the case K < and L < the partition function of this Hamiltonian can then be 
transformed to the partition function of an effective Ginzburg-Landau-Wilson-Hamiltonian for which one can show 
that it has the critical behavior of the 3d-XY model. 

A conne ctio n with experimental results obtained for actual crystals is established by applying the method described 
in section VT1 to the DIS model. For BCCD, for example, this leads to a descriptiono which is at least as good as 
the description by the ANNNI model, the latter being a special case of the formeiO. 



D. Interrelations of various pseudo spin models 

A general r-spin model containing r spins per cell cr*, . . . , er£ can also be described as a q = 2 r -state model. The; 
state variable p n = 1, . . . , 2 r and the spin variables g_ n = (cr* , . . . , c^) are, for example, connected by the equationEil 



2 

i=l 



Special cases of r-spin models are, e.g., 2 r -state Potts models. 

Pseudo spin systems in the narrower sense are described by Hamiltonians of the form 



U = V A ij cr* a 3 ' 



For uniaxially modulated systems with nearest neighbor interactions in the planes perpendicular to the direction of 
the modulations this Hamiltonian may be rewritten as 



<n,m> 
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The matrix V denotes couplings perpendicular to the modulation direction, the sum running over all nearest neighbor 
pairs (n, m) in the planes, whereas W denotes couplings along the direction of the modulations. 



In the previous sectio ns we discussed three different pse udo sp in models: the ANNNI model (section VIA), the 



AANNDI model (section VI B )) and the DIS model (section VIC). Whereas the ANNNI model is an one-spin model, 
the DIS and the AANNDI model are, in their original formulation, two-spin models. However, in the treatment of 
the AANNDI model Kurzyhski considered three independent order parameters (t a ), (& A ) and (^ A ) — (t a <t a ), i.e. 
the AANNDI model is treated as a model with the three pseudo spins t a , <j a and ui A . The DIS model is, on the 
contrary, treated as a model with two pseudo spins r and a per cell. 



Both the AANNDI and the DIS model can be derived from the general Hamiltonian (12). 

The DIS model represents a model which is somewhat more general than the ANNNI model but nevertheless so 
simple, that it allows the derivation of explicit results (dependent on a few well defined local interactions) which can 
be compared to experimental data. This powerful model contains some other much discussed models as special cases: 

Since a general r-Ishaa spin model (i.e. a model, where r spins describe the state of every unit) can also be considered 
to be a 2 r -state modelEil (i.e. a variable p = 1, . . . , 2 r describes the state of every unit), it follows that the DIS model 
can be formulated as a four-state model. For the special case K = L < and J = J' < the DIS Hamiltonian can 
be rewritten in the forrrO 



n 



-Co ^2 COS (Pn - Pm)) - C C ° S (f ( Pn ~ Pm + A ) 



V 1J 

<n,m> <n,m> 



where we introduced the variables 



Co := -2J , C := VM 2 + 4K 2 and A := -- arctan ( — ] . 

7T V 2K J 

This is exactly the Hamiltonian of the four-state chiral clock (CC4) modeS. It follows from the exact expansion of 
the free energy for the special case K = L < of the DIS model that in the limit, T — > the CC4 model exhibits 
a complete devil's staircase. This is in marked contrast to an earlier statement^ that only specific phases spring 
from the MP. A detailed analysis of this rather complicated calculation revealsEj that in the earlier investigation the 
treatment of the in-layer interactions was erroneous, thus leading to wrong expressions for the free energies. In fad- 
recent Monte Carlo simulationsE3 of the CC4 model show the existence of modulated phases not predicted by ref.E£l 
but in accordance with oup_results. 

It can further be showrO that a special case of the DIS model (two pseudo spins per cell with nearest neighbor 
interactions) can be mapped exactly onto the ANNNI model (one pseudo spin per cell with nearest and next nearest 
neighbor interactions). At the first sight this is surprising, since with respect to the range of direct interactions the 
ANNNI model is more general. For the special case considered no direct couplings between neighboring cr-spins are 
retained, i.e. L = and J' = in cqn. (pL4[) . The cr-spins then mediate an indirect next nearest neighbor coupling 
between the r-spins. The resulting relation between the coupling parameters of the two models is temperature- 
dependent. Its explicit analytic form istil 

Jo = — J, 

h = -K, 

J 2 = In cosh GSM). 
zp 



VII. CONNECTION OF MODEL CALCULATIONS WITH EXPERIMENTALLY DETERMINED 

PROPERTIES OF ACTUAL MATERIALS 

A. Theoretical derivation of phase diagrams in terms of experimental control parameters 

All the reviewed theoretical models for the description of uniaxially modulated materials have in common that the 
displayed phase diagrams are twodimensional sections of a higher dimensional space spanned by the model parameters 
(for example reduced temperature, reduced interactions etc.). The experimental phase diagrams, however, are spanned 
by external quantities like temperature, pressure or applied fields. It is thus necessary to find a mapping between 
the experimental quantities and the theoretical parameters in order to compare the experimental and theoretical 
phase diagrams. In the following we will present a procedure for transforming a two-dimensional theoretical phase 
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diagram from model parameters, la temperature-pressure variables. This mapping has originally been formulated for 
the ANNNI and the DIS modelfl'El in order to describe BCCD. Since the experimental wave numbers q are given as 
multiples of 2ir/c, the assignment of units to half cells yields (for the ANNNI model) q — vj (Z\ + Z2 + ■ ■ ■ + Z v ) for the 
relation between the phase symbol (Z1Z2 ■ ■ ■ Z v ) and q. We will now present a more general approach transforming 
arbitrary theoretical phase diagrams into temperature-pressure phase diagrams. 

The method requires the knowledge of the thermal expansion coefficients, the elastic constants, and information on 
the critical line T cr ;t (p C rit ) separating the normal phase from the lower-symmetry phases in the T-p diagram. From 
the determination of the critical wavenumber for different pressureso a fit <? C rit(]?crit) can be found. 

Let the two quantities spanning the considered two-dimensional theoretical diagram be k and 9. In the ANNNI 
model, for example, k is given by k — — ^f- whereas 9 is the reduced temperature 9 = ^fp- The critical line # C rit( K ) 
and the critical wavenumber gfcrit(«) result from an analytical calculation (ANNNI, AANNDI, DIS model) or from a 
numerical fit (DIFFOUR, Chen and Walker's model). 

In the first step the critical line #crit(K C rit) ^ s transformed into a critical T-p-line. For a given k = kq the critical 
values gcrit(Ko) and C rit(Ko) are obtained. On the other hand, as the wavenumber q cr n is now known, the values 
Pcrit(<7crit) and T cr it(y>crit(9crit)) for the critical pressure can be derived, thus leading to a theoretical T-p critical line. 

In the next step we consider an arbitrary point (n,9). The temperature coordinate is obtained by the simple 
transformation rule ksT = fc ^ Tcrit 0. Depending on the models this rule may be replaced by a more sophisticated 
approach. For simplicity, we assume that thq-.pijodel interactions depend only via the mean lattice constant on 
temperature and pressure, i.e. Jk = Jk (o(T , ,p))Qc3. Under this assumption the interactions are constant along lines 
where a is constant. Lines with constant a are especially lines with constant volume. Expanding the volume about 
its value for a given temperature Tq and pressure po leads to the equation 

dV dV 
V(T,p) = V(T a , Pa ) + -^(T , Po ) ■ (T - T ) + ^(To,Po) ■ (p - Po) 

ld 2 V 2 1<9 2 U 2 

+ 2Qf2( T °'P^ ■ (T-Tv) + 2 9^( T »>W) ' (P-Po) 

d 2 V 

+ dfdp <yTo,Po ^ ' ( T-T °) (P~P°) 

or 

= V(T,p)-V(T ,po) 
V(T ,p ) 

= 3a ■ (T - T ) - k ■ (p - po) + • (T - T f 

~ 2 d~ P ' {p Pa) -of^ T - To ^p- ( 15 ) 

where we introduced the thermal expansion coefficient a — and the compressibility k — —y^- The value 

of v is determined by the intersection point of the line of constant volume with the para phase boundary, i.e. v = 
v(r erit ,porit)-V"(T ,po) ^ i nser t m g ^ ne temperature T in eqn. fll5| ) finally leads to the respective pressure. 

Figure shows the theoretical temperature-pressure phase diagram derived from the ANNNI model for the de- 
scription of BCCDEa. This diagram shows the same topology as the experimental diagram, similar extensions of the 
modulated phases, and the observed branchings. The phase diagram for the DIS model is transformed in a similar 
way as the ANNNI-phase diagramE3. 

The observecD increase of transition temperatures with uniaxial compressional strains —en corroborates the as- 
sumption that the interactions depend on p and T only via the lattice constants. 

This approach especially allows the determination of the effective couplings-«A. as a function of the lattice parameters. 
They can also be derived from the shape u>(q) of a softening phonon brancl£3. A linear chain with up to Ath nearest 
neighbor couplings with dispersion relation 

N 

[hio(q)f = ^ J k (1 - cos kq) 
fc=l 

is considered as analogon to the crystal. A least-squares fit of this relation to the soft phonon branch for various 
temperatures T above the transition to the modulated phases yields the dependence Jfe(T). 
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B. Structure predictions in terms of the SALM model and its derivatives 



For uniaxial modulations, the mean values of atomic displacements do not vary in planes perpendicular to the axis 
of modulation, which is c in the case of BCCD. Hence, the generalized variables have (at the thermal average) values 
Ql = Q e (3), Ql = Qe+i(3) and Qt = Qt{4), Qi = <3f+i(4) where £ = 2n 3 and £ + 1 = 2n 3 + 1 label consecutive 
layers (dotted vertical lines in figure g) formed from quadruples of blocks belonging to cells at positions (711712713). The 
displacements related to the , layer variables' Qe(3) and Qe(A) are so-called , layer modes'. The notation is adapted 
to the fact that these layers can be considered equivalent as the space group operation Q 5 — {C2I5I3} generates a 
pseudoperiodicity of half a lattice constant along the c-axis by transfo rmin g one layer into the next. This construction 



yields , layer variables' as in Chen and Walker's model (cf. subsection |VB| ). Specific differences arise for modulations 
along c (e.g. BCCD): The layers in Chen and Walker's model are located at positions zg = (2£ + j) and Z( + i = (2£+ 1) 
and thus intersect (in our nomenclature) blocks belonging to different quadruples, as can be seen from figure [l] and 
figure U (solid vertical lines). Chen and Walker's choice of layer modes hence is a linear combination of the layer 
modes of the SALM model for modulations along c. Two consecutive nitrogen atoms (along the c-axis) are displaced 
simultaneously by each of the two modes for the respective layer they belong to. For both models, the first layer mode 
displaces both atoms in the same direction along b whereas the ^-displacements produced by the second layer mode 
are in opposite direction. Depending on how one combines the nitrogen atoms and the layers, either both modes are 
necessary for a satisfactory approximation of the displacements (Chen and Walker's model) or only one (the layer 
mode related to Qi(3) in our nomenclature). Hence, we tried a structure prediction for the displacements of the 
nitrogen atoms in the fourfold phase of BCCD either in terms of the layer modes of Chen and Walker's model or 
in terms of relevant cell modes. If, following Chen and Walker, two consecutive nitrogen atoms are assigned to the 
planes given by full vertical lines in figure ||, their displacements are quite different and two modes are needed for a 
reasonable description. If, however, the nitrogen atoms are assigned to the planes given by dotted lines, as proposed in 
the present paper, their displacements are very similar and one (relevant) mode gives already a very good description 
as shown in figure [| An application of the DIS model to BCCD, which corresponds to taking into account a second 
set of SALMs as opposed to the ANNNI model with one set of SALMs, yields even better agreement of calculated 
and measured structures. 

If both the ANNNI model and the DIS model are treated in mean field approximation, there is no need in dis- 
tinguishing pseudo spins derived from amplitudes of cell modes for front or back half cells: in an application of the 
ANNNI model to BCCD the pseudo spin a is derived from Q 3 ,Q 7 . In the DIS model the pseudo spin r is derived 
from Q 3 , Q 7 whereas the pseudo spin a is derived from Q 4 ,Q 8 . 

Since the pseudo spins a (for the ANNNI model) or a and t (for the DIS model) stand for the amplitudes of certain 
local modes, their transformation behavior under space group operations can be determined and relations between 
the couplings can be found. It can be shownEJ that the antisymmetric coupling M between the two pseudo spin 
subsystems of the DIS model is not introduced ad hoc but follows necessarily from symmetry. 

Making use of the transformation behavior of the generalized variables, expressions for the spontaneous polarization 
of the crystal in terms of the amplitudes Qe(3) and Qi(4) can be derived. Since the components of the spontaneous 
polarization vector Pg must transform like any vectorial quantity under the operations of the space group Prima, 
only such terms are allowable in an expansion of P5 in terms of relevant mode amplitudes that are in accordance with 
symmetry. Expanding P5 up to second order in Qe(3) and Qe(4) and keeping only terms involving same or adjacent 
layers, the following expressions can be derived: 

Ps,x = ^(-l)V^(3)^ + i(3) + PtQM)Qt+i(±) + PiQe(3)Q e (A) 
(Q/(3)Q /+ i(4) - Q e (A)Q e+1 (3))] + ... 

Ps, v = J2(- 1 ) ip vW) + --- > 

1 

Ps,z = ]T P^Qt(3)Qe(4) + ... . 

£ 

P*, ■ ■ ■ , P£, Py and Pi are constants specific for BCCD. A sequence of equal Qi leads to alternating signs of x- and 
^-displacements in consecutive layers as the respective layer modes are transformed one into another by the screw 
axis G 5 - This is as a consequence of the proposed definition of the local modes. Thus, in the case of BCCD, a 
correspondence of pseudo spin averages Sg to Qg(3) with alternating signs was performed in our previous worko. The 
corresponding expressions for the DIS model with 

Q0) — > {-!)%, 
Q t {4) — > (-I)'** 
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read 



Ps,x = ^(-l) £ [PiiA+l + Pjs^+i + P^Siti + (t t B t+1 - Site+l)}, 

e 
i 

P S ,z = J2 p * ttst - 



(16) 



For the values of pressure and temperature applied in measurements in bcceB 
izations are obtained depending on the wavenumber: ^ = : Ps,x\ v ~ ' xt " 



the following spontaneous polar- 
Ps, v ; ~ = no spontaneous 



even - - 1 a > Xi n odd 

polarization. This is in perfect agreement with the spontaneous polarizations obtained from the DIS model for the 
corresponding model parameters and the assumption that a symmetry mode of A3-symmctry gives the primary con- 
tribution to the structural deformation. If only one mode per half cell is retained (Q 3 , Q 7 ), the expressions for Ps,x, 
Ps, y and Ps, z are gwen by the terms in eqn. ( |l6|) depending only on t(, in complete accordance with the terms 
we proposed in ref.E3. Recently new dielectric investigations on BCCD revealed the so-called Ts-anomalytZl in the 
regions of the commensurate \- and ^-phases exihibiting the characteristics of a commensurate-commensurate phase 



transition where the wave number is unchanged 
model; a detailed discussion will be given elsewhcrel 



same phenomenon occurs in structures derived from the DIS 



VIII. THE INFLUENCE OF POINT DEFECTS 



In subsection VII A it was shown that - even if there is not yet an ab initio type of theory - it is possible to 
establish a quantitative connection between experimentally determined T, p-phase diagrams and theoretical diagrams 
formulated originally in terms of model parameters. In the case of an ANNNI model description, the latter are the 
couplings Ji, defined, as explained above, by thermal averaging over local modes. A simplifying step in this procedure 
was the plausible assumption that these model parameters depend on T and p mainly via the total strains or via the 
lattice constants. By making use of macroscopic material properties like thermal expansion, elastic coefficients etc., 
this lattice constant dependence of the couplings Jj was explicitely derived. The similarity of the theoretical and the 
experimental phase diagrams support this procedure. The success of this approach made it promising to interpret 
defect influences, e.g. on phase diagrams, along the same lines. For the general treatment we shall extend the basic 
assumption somewhat, allowing for a dependence of the effective interactions not only on T and dilatation but on 
T and general strain components. We shall first formulate this extended version, then include defect contributions 
to the strains, and finally give some explicit results on defect influences. The method is exemplified for the ANNNI 
model; it can, however, easily be transposed to other models. 

Let the phase diagram be formulated in the space of the experimental control parameters temperature T and 
components a E of the tensor <r E describing a homogenous external (applied) stress field. The aim of the method 
to be described is to obtain an approximate theoretical statement as to which one of the many modulated phases 
occuring in the model will be stable at a point (T,a E ) in parameter space. 

If results of ab initio electron theory or of atomistic interaction potential theory followed by a thermal averaging 
procedure were at hand, these would yield the effective interactions Ji(T, <r E ) and the stable phase could directly be 
read (T,a E ) from the theoretical ANNNI phase diagram in (T, Jo, Ji, J2)-space for any set of control parameters. 

If the dependence Ji(T, a E ) is not known from basic theory, the following procedure, which makes use of macroscopic 
properties and characteristics of the para phase boundary (PPB) yields approximate results: let the empirically 
determined PPB and the modulation vectors on the PPB be given by T = T(cr E ) and q = q{cr E ) and assume the 
effective model parameters Ji to depend on T and cr E (only) via certain linear combinations of the stress tensor 
components Eif 

Ji Ji (Pynn^'* n 7i) — • JiiS^ 1 

(summation over repeated subscripts). The coefficients Rmn are chosen in a plausible way taking the lattice structure 
into account. The displacements are counted from the para phase positions at some temperature To above but close 
to the PPB, zero applied stresses cr E and no defects. For crystals free of defects the total strains are, in linear 
approximation, the sum of elastic and thermal expansion contributions: 

£mn = S mn kjO- E j + a mn (T — Tq) (17) 
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and 



Q« = R$ n S mnkj a% + Rgla mn (T - T ) =: Ajgag + B^(T - T ). (18) 

For this step, values of the macroscopic elastic moduli and thermal expansion coefficients are required. Relation Jl8| ) 
may be (and was) improved by taking higher order terms and coefficients into account. For any point (defined e.g. by 
the external stress ar E ) on the PPB, the corresponding set of model parameters Ji may be determined by observing 
that in the ANNNI model the modulation vector on the PPB depends only on the interaction ratio k = —J2/J1 [cf. 
eqn. @] : 

ffcrit(t) = tt - arccos-^- for k > 0.25. 
2tt 4k 

Thus k is known for any point on the PPB. Now the model PPB is given by 

kT clit = 4J + (2k + -L)Ji =: 4J + /(k)Jl 
4k 

Jo and J\ can be determined from the (experimentally determined) dependence of T cr - lt on /(k) for variations of <t e 
which leave 

fi« = Ajgo$ + B«(T crit (<x £ ) - T ), i = 0,1 

and thus also Jo and J\ invariant. Once k, Jo, Ji are determined, J2 is obtained from k and Ji. In this way one 
obtains corresponding values of Ji and that is a numerical representation of the function Jj(fiW). One can 
thus assign to any point (T, cr E ) (also to points off the PPB) the corresponding values Ji. Reference to the ANNNI 
(T, Jo, Ji, J2)-phase diagram then yields the phase stable at (T, cr E ). 

For defect crystals the contributions e D of the if"j)int defects to the local strains have to be included in eqn. (|T^). 
The deformation due to a point defect (see e.g. refOO) is described by its force dipole tensor Pij. With the material's 
clastic Green's function Gjy(r, r'), which is expected to be, in good approximation, the same for different modulated 
phases, the strain field at r produced by (equal) defects at positions v" is 

v 

and that produced by an uncorrelated distribution of defects with local defect density p D {r) is 

e° n (r,p D )= J G m ^ n (r,r')P kj dr' (19) 

with f t k' ■— df/dxk'- In extension of the procedure described above, the local model parameters, e.g. the ANNNI 
couplings Ji (r) are then assumed to depend on the local strains e and via thesejj] not only on temperature and external 
stresses cr E (e.g. on pressure), but by eqn. (|i~9|) also on the defect distribution p D (r): 



J(r) = J (e(r)) = J ({S mnkj a^ + a mn (T - T ) + e£„(r, P D )}) 
Dividing the defect strain into its average and a position dependent contribution one has 

e D (r) =^V)+ (e D (r,p D )-7^Y 



The influence of the average strain corresponds to an extra term aj^ = Ckj pq £p q , which has to be added to the external 
strain a E j . From the experimentally observed average deformation e connected with a defect concentration p D , the 
dipole tensor can be determined. One has 



In addition, the interactions (between two units) will be changed quite apart from the defect induced strain simply due to 
the fact that in one of the two units there is a substitutional or an interstitial atom. This effect is not considered here. 
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If our basic assumption is correct that the model parameters Jt depend on temperature and external stresses mainly via 
the total strain field including elastic and thermal contributions and if the main defect contribution is a function of the 
average defect induced strain e D (p D ), then the dominant impurity effect on the phase diagram should be an alteration 



obtained by replacing the external stresses a^- by <jjy + p D Pkj = a^j + ^kjmn^mn- This relation was checiq 



IOC 
10l|-103 



for 
for 



brominated BCCD under hydrostatic pressure for which phase diagrams were experimentally determined 
various Br-concentrations. Already the rough approximation of the variables £l(i) by the-dilatation e mn /3, i.e. the 
average relative volume change, leads (with experimentally determined material constantse2-3't2j) to a good agreement 
with experimental data: For Br-concentrations of 1.8%; 7.3%; 20% the calculated pressure shifts are -20 MPa; -81 
MPa; -222 MPa as compared to experimentally observed shifts of -18 MPa; -R1 MPa; -238 MPa respectively. 



The influence of spatial variations of the defect strains was investigated ! 100 ! for a situation with simple model char- 



acteristics, control parameters and defect arrangements: the coefficients R in fiW were chosen to define the dilatation 
part of the strains {Rmn = 5 mn ), the external stress is a hydrostatic pressure (<r,y = —pSij), and the defect distribution 
is such as to produce a dilatation strain field e D := ej? which varies only in one direction (perpendicular to the mod- 
ulation axis), with a position (x) dependence described by a saw tooth function. The nearest neighbor interactions 
Jo and Ji were set equal to each other. Depending on the amplitude of the spatial variation of defect concentration 
or strain, this scenario leads to a coexistence of several phases for given values of temperature and pressure. A 
small section of the phase diagram is considered and the general procedure described above for calculating the phase 
diagram is modified correspondingly: the strain field and the free energy densities p of the phases 7 are expanded 
about values of temperature, pressure and position corresponding now to a point in the middle of that section at a 
local position in the middle of a saw tooth: 

P{J i [n(T,P,e D (x))lT} = ft + AT, AT + A^Ap + AJAx 

with / 7 = p{T Q ,p Q ,x°) and 

_ dpdJidn dp ^dp^dJidn _ dp dj, an de D 

T ' dJ, dVLdT dT' p ' dJi dn dp' x ' dJi dn de D dx ' 

With the position x 1 of the boundary between phases 7 — 1 and 7, and with 7 = the first and 7 = m the last phase 
to occur, the total free energy is 

F(p,T,{x^}) = J2 Pdx. (20) 

Minimizing F with respect to the boundary positions x 1 yields the local extensions of the different phases as function 
of T and p. . — . 

Explicit calculations!^ were performed for the section of the phase diagram show in figure ^. The temperature 
(or pressure) dependence of physical quantities which are different in different phases (i.e. which depend on the 
modulation wave number q) will vary in a typical way with the defect concentration gradient. As a representative of 
such quantities we consider q itself. Its average q is obtained from the values x 1 , and the wave numbers q 1 . 

Results of such calculations are shown in figure |l^. For crystals with a spatially constant defect density, q(T) 
is characterized by a series of plateaus corresponding to one modulated phase each with discontinuous transitions 
between them (figure |Io|a). Upon introducing a spatial (saw-tooth-like) variation of the defect density, the plateaus 
become gradually narrower and the transitions become continuous. This is shown in figure G5 the amplitude of 
the defect density variation increases in equal steps from figure ^fj| a (amplitude zero) to figure px|f. The transitions 
intervals now correspond to the coexistence of two or more phases whose boundaries are continuously shifted with 
temperature. For sufficiently high concentration gradients some plateaus disappear completely, indicating that the 
stability range of the corresponding phase 7 is so small that there is no temperature for which 7 is the only stable 
phase in the whole specimen. 



IX. DISCUSSION AND OUTLOOK 



The theoretical description of materials exhibiting modulated phases may be subdivided into four main points: 
1) Eventually the description should be based on an ab initio electron theory (e.g. rei!23) or, if that is too difficult 
to handle with the necessary accuracy-for the rather large number of atoms per unit cell, on calculations using 
empirical atomic interactions (e.g. reill23). A suitable procedure would be to determine total energies for the para 
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phase configuration and for a number of appropriately chosen deformations. The respective results should allow the 
selection of adequate simplified models, the formulation of the respective effective Hamiltonians, and a determination 
of the model parameters. These models should then be handled as described under points 2), 3), 4). So far there are 
scarcely any results of this kind. Calculations along these lines, however, would allow to go, for a not too complicated 
system, all the way from an atomistic theory to explicit statements on phase diagrams, polarizations etc. In such a 
treatment approximations are of course necessary and should be made (in the original atomistic theory as well as in 
the consecutive model calculations); but one should avoid ad hoc assumptions (fits) as far as possible. 

2) The second part of the description is the formulation and investigation of the properties of model systems with 
one or a few relevant local variables like e.g. the DIFFOUR model. Especially if this is the first step (i.e. if one cannot 
deduce these models explicitely from one of the ab initio theories of point 1), all symmetry information should be 
taken into consideration in the formulation of the models. The main task is the statistical mechanics treatment which 
should start from the respective (effective) Hamiltonian and should yield information on phase diagrams, orders 
of the transitions, polarizations etc. The calculations may be and have been performed either in a more or less 
analytic way by using schemes of approximation like e.g. the mean field or the self consistent phonon method. Or, 
alternatively, numerical approaches (molecular dynamics simulations) were used. So far, calculations for such models 
were performed mainly for very simplified geometries like linear chains. 

3) A further step in the simplification process is the projection of models with a (small) number of continuous local 
variables onto models with discrete states, or, finally, onto pseudo spin models. Depending on the number of relevant 
original local modes one obtains pseudo spins with one or more components per crystallographic unit. Substantial 
advantages of this type of description are, that the powerful apparatus of Monte Carlo simulations can be applied and 
that analytic approximations exist which allow the explicit calculation of many interesting physical properties. This 
part of the theory is well developed. This holds especially for mean field treatments which in some instances where 
compared with more pretentious methods and shown to yield a good description. 

4) Finally there is the problem of translating the results of model calculations obtained as functions of the model 
parameters into data which can be quantitatively compared to experimental data, i.e. which are expressed as functions 



of experimentally given quantities like pressure or stresses (besides temperature) . In section VII it was shown how this 
can be done. By making use of macroscopic material constants, a physically motivated mapping from one parameter 
set to the other can be performed. For the material considered (BCCD) the results show a satisfactory agreement 
with experimental data. 

In the last years considerable effort has been addressed to the approaches described in section VI, that is to the 
formulation of various p-state and pseudo spin models and to their statistical mechanics. Methods applied were MFA 
calculations, exact series expansions for selected regions in parameter space, and Monte Carlo simulations. Combining 
these approaches, a reliable theoretical determination of most of the interesting thermodynamic properties of such 
models is possible. 

A transformation of these model proper ties in to a form directly comparable to experimental data can be performed 
by the procedure described in subsection VII A or by a similar scheme. An interesting by-product of this method is 
information on the dependence of the model parameters on the experimentally given state variables and especially 
on the lattice distortions. Combined with a consideration of experimentally determined atomic positions in some of 
the more stable phases, this should further the understanding of the local dynamics responsible for the transition 
sequences. 

Of special interest, finally, would be the determination of model parameters - for not too complicated materials - 
from an atomistic theory. It should be possible, with the nowadays available computing power, to carry through the 
necessary calculations. Investigations of this type would clarify the details of the relevant atomistic mechanisms and 
e.g. help to understand why the phase diagrams of members of the same homologous family of materials may differ 
drastically. 
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FIG. 1. The blocks i = 1 ... 8 belonging to the unit cell (depicted as wireframe) at the origin O. 



FIG. 2. Mean field phase diagram of the ANNNI model with J = Ji. 

FIG. 3. The different fixed point types. Shown are the positions of the four eigenvalues of the matrix G k in the complex 
plane. 

FIG. 4. Ground state of the DIS model. The pair of upper/lower signs per column of the symbols represent the signs of the 
r- and cr-spin averages in consecutive layers: I is the fully ferroelectric phase, II the fully antiferroelectric phase, whereas III 
and IV are mixed phases in which one spin type orders ferroelectrically, the other antiferroelectrically. The structure of the 
phase V repeats itself after 4 layers. The vertices are given by (A, |y) = (1,0), (0,1), (-1,0) and (0,-1). 

FIG. 5. Global phase diagram of the DIS model with k_ = —0.1 and ji = jj = —0.25. Only few modulated phases 
are shown. 6 = ^jj- is the reduced temperature, k± = K ^j L , and K, L, M are the different couplings in direction of the 
modulations. 

FIG. 6. Global phase diagram of the DIS model with A = £, = 0.1 and jr = C = -0.25. 



FIG. 7. (p, T)-phase diagram for BCCD calculated from the ANNNI model. Dashed lines: boundaries to higher commen- 
surate or incommensurate phases. 

FIG. 8. Measured (o) and calculated (•) ^-displacements of the nitrogen atoms in the fourfold phase of BCCD. The mod- 
ulation profile was calculated from the ANNNI model, i.e. with one variable per unit (crystallographic half cell). Solid and 
dotted lines: centers of layers in Chen and Walker's model and in the SALM model respectively; consecutive nitrogen atoms 
are assigned to the layer between them. 

FIG. 9. Section of the ANNNI model phase diagram exhibiting some of the more stable phases for which the results shown 
in figure [l^ were obtained. This section is part of the region occupied by the phases observed in BCCD. 

FIG. 10. Averaged modulation wave number q vs. O = kT/Ji. Plateaus correspond to phases (56) (1), ((56) 2 556) (2), 
(56556) (3), (56(556) 2 ) (4) and (556) (5). p D x and a tX are increasing from a) (a = const) to f) in equal steps. In all six cases 
the sum of applied pressure and chemical shift at Xq is kept constant (e.g. by keeping p D (xo) and a(xo) fixed). 
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